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ABSTRACT 

The motion and variability of the radio components in the low mass X-ray binary system 
Sco X-l have been monitored with extensive VLBI imaging at 1.7 and 5.0 GHz over four years, 
including a 56-hour continuous VLBI observation in 1999 June. We detect one strong and one 
weak compact radio component, moving in opposite directions from the radio core. Their relative 
motion and flux densities are consistent with relativistic effects, from which we derive an average 
component speed of v/c=0.45±0.03 at an angle of 44° ±6° to the line of sight. This inclination of 
the binary orbit suggests a mass of the secondary star that is < 0.9 M@, assuming a neutron star 
mass of 1.4M Q . We suggest that the two moving radio components consist of ultra-relativistic 
plasma that is produced at a working surface where the energy in dual-opposing beams disrupt. 
The radio lobe advance velocity is constant over many hours, but differs among lobe-pairs: 0.32c, 
0.46c, 0.48c, and 0.57c. A lobe-pair lifetime is less than two days, with a new pair formed near 
the core within a day. The lobe flux has flux density that is variable over a time-scale of one hour, 
a measured minimum size of 1 mas (4 x 10 8 km), and is extended perpendicular to its motion. 
This time-scale and size are consistent with an electron radiative lifetime of < 1 hr. Such a short 
lifetime can be caused by synchrotron losses if the lobe magnetic field is 300 G or by adiabatic 
expansion of the electrons as soon as they are produced at the working surface. The lobes also 
show periods of slow expansion and a steepening radio spectrum. Two of the core flares are 
correlated with the lobe flares under the assumption that the flares are produced by a coherent 
energy burst traveling down the beams with a speed > 0.95c. 

The radio morphology for Sco X-l differs from most other Galactic jet sources. Possible 
reasons for the morphology difference are: Sco X-l is associated with a neutron star, it is 
a persistent X-ray source, the source is viewed significantly away from the angle of motion. 
However, the lobes in Sco X-l are similar to the hot-spots found in many extragalactic radio 
double sources. Scaling the phenomena observed in Sco X-l to extragalactic sources implies 
radio source hot-spot variability time-scales of 10 4 yr and hot-spot lifetimes of 10 5 yr. 

Subject headings: binaries:close; galaxies: jets; radiocontinuum:stars; stars:neutron; stars:individual(Sco 
X-l); X-rays:stars 
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1. Introduction 

Since Sco X-l was first detected at radio fre- 
quencies in 1969 (Abies 1969), many observations 
at radio, optical and X-ray frequencies have been 
made in order to understand the physical pro- 
cesses associated this object. It is identified with 
a 13-mag binary system with an orbital period of 
0.787d (Gottlieb et al. 1975). The degenerate ob- 
ject is probably a neutron star, and the companion 
is an unknown spectral type with about a solar 
mass (Cowley & Crampton 1975). Sco X-l is one 
of the most intense persistent X-ray sources and 
one of about 30 presently known low-mass X-ray 
binary (LMXB) systems. It is a 'Z-type' LMXB, 
named for a characteristic shape of its X-ray 
color-color diagram. It also exhibits quasi- 
periodic X-ray oscillations (QPO) from a few 
Hertz to kilo-Hcrz (Hasinger & van der Klis 1989; 
van der Klis et al. 1996). 

Radio observations in 1970 detected two sym- 
metrical radio sources about one arcmin from 
Sco X-l (Hjellming & Wade 1971). They were 
thought to be radio lobes of Sco X-l, but 
were later shown to be stationary extragalac- 
tic background objects having no relationship 
to Sco X-l which has a significant measured 
proper motion (Fomalont & Geldzahler 1991). 
The results from a multiwavelength campaign 
(Hjellming ct al. 1990) also suggested that there 
was only weak correlation of the radio emission 
with the optical and X-ray emission. 

High-resolution radio observations at 5.0 GHz 
with the Very Long Baseline Array (VLBA) 
were initiated in 1995 in order to determine 
the trigonometric parallax of Sco X-l. From 
eight VLBA observations through 1998 August, 
a distance of 2.8 ± 0.3 kpc was determined 
(Bradshaw et al. 1999). This distance (previ- 
ous estimates ranged from 200 pc to 2000 pc 
(Westphal et al. 1968)) established that the X- 
ray luminosity was near the Eddington limit 
as inferred from X-ray models of Sco X-l and 
other z-type LMXB's (Lamb 1989; Penninx 1989; 
Vrtilck et al. 1991). The radio emission during 
these VLBA observations showed variability in in- 
tensity and structure over periods less than one 
hour, with discrete components moving away from 
the core at relativistic velocities. 

With such variable radio emission in Sco X- 



1, this source was established as another ex- 
ample of a Galactic-jet variable radio source 
(Mirabel & Rodriguez 1999) with strong simi- 
larities to the jet and lobe phenomena associ- 
ated with luminous radio galaxies and quasars 
(Blandford & Rees 1974). The evolutionary time 
scales associated with accretion disks and jets 
scale roughly with the luminosity and mass of 
the degenerate object. Thus, significant changes 
should occur over hours and days in the evolution 
of Galactic sources that would take millions of 
years in quasars. 

Experience gained from the eight VLBA obser- 
vations used to determine the parallax of Sco X- 
1 showed that the four to six hour observations 
made every six months were not sufficient to deter- 
mine a coherent story for the radio emission from 
Sco X-l. Each observation was too short to follow 
much of a component's evolution, and the indi- 
vidual observations were too separated in time to 
follow the evolution. For this reason, we observed 
Sco X-l with VLBI resolution for a continuous 56- 
hour period in 1999 June. The results from this 
series of observations were crucial to understand- 
ing the radio evolution of Sco X-l and formed a 
framework for understanding the previous VLBA 
observations. 

We have already reported results associated 
with some of these data. The parallax was 
given by (Bradshaw et al. 1999), and Paper I 
(Fomalont et al. 2001) discussed the kinematic 
properties of the radio source; the space motion 
of the lobes and the speed of the energy flux in 
the beam. This paper is organized as follows. The 
nine VLBI observations, and the data reduction 
methods are described in §2. Because we display 
VLBI snapshots (images made from observations 
less than about two hours), the difficulties and 
ambiguities with such imaging are described. The 
images in §3 from the 1999 June 11-13, 56-hour 
VLBI observation and the 1998 February 27-28, 
observations are presented and discussed. Then, 
images from the other less extensive observations 
are shown and compared with the more extensive 
observations. In §4 we review the kinematic prop- 
erties of Sco X-l (Paper I) and estimate the stellar 
component masses of Sco X-l from the orbital in- 
clination. The detailed physical properties and 
suggested production mechanisms associated with 
the lobes are discussed in §5. The optical, X-ray 
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and radio properties and evolution of the core will 
be discussed elsewhere (Geldzahler et al. 2002). 
The comparison with extragalactic radio sources 
is given in §6, and our conclusions are summarized 
in §7. 

2. Observations and Data Processing 

Table 1 contains the log of all VLBI observa- 
tions of Sco X-l. The first eight observations 
with the VLBA, between 1995 August and 1998 
August, were planned to optimize the determi- 
nation of the trigonometric parallax. The 56- 
hour monitoring of Sco X-l in 1999 June was 
composed of a series of seven consecutive eight- 
hour observations amongst three different VLBI 
arrays: VLBA+VLA (the Very Large Array), the 
APT (Asia-Pacific Telescope), and the EVN (Eu- 
ropean VLBI Network + additional telescopes). 
Extensive optical and X-ray observations were also 
made during this period (Titarchuk et al. 2001; 
Bradshaw et al. 2002; Geldzahler et al. 2002). 

2.1. Observations and Data Correlation 

The VLBA observations alternated between 
Sco X-l for 9 minutes and the calibrator (1504- 
166) for 2 minutes. Three pairs of observations 
were made at 5.0 GHz followed by one pair at 
1.7 GHz. The APT and EVN observations were 
made solely at 5.0 GHz and the observing sched- 
ule alternated between Sco X-l for 19 minutes and 
the calibrator source for 3 minutes. The VLBA 
and EVN observations were observed with a to- 
tal bandwidth of 128 MHz with one-bit sampling, 
whereas the APT observations used a total band- 
width of 32 MHz, but with two-bit sampling. The 
VLA was used in two ways: as one element of the 
VLBA and as a stand-alone array from which arc- 
second resolution images were obtained for Sco X- 
1 (Geldzahler et al. 1983; Bradshaw et al. 1997). 

A background radio source, about 70" NE of 
Sco X-l, has a correlated flux density of 8 mJy at 
5.0 GHz and 14 mJy at 1.7 GHz. This object is 
the north-eastern of the two sources near Sco X- 
1 which were previously thought to be related to 
the binary system (Fomalont & Geldzahler 1991). 
This source is sufficiently strong and compact to 
be used as the primary phase calibrator for Sco X- 
1 (Bradshaw et al. 1999) at both frequencies. We 
will designate this radio source as the North-East 



in-beam calibrator (NEIBC). Because the effective 
primary beam of the phased arrays (VLA, West- 
erbork and ATCA) were smaller than the 70" sep- 
aration between Sco X-l and the NEIBC, these 
arrays alternated observations between Sco X-l 
and the NEIBC with three-minute cycles. All 
other telescopes pointed midway between the two 
sources and observed them simultaneously. 

The VLBA and EVN data were processed by 
the VLBA correlator in Socorro, NM, USA. The 
data from the APT observations were processed 
with the S2-correlator in Penticton, BC, Canada. 
All data associated with the observations of Sco 
X-l were correlated at two positions: (1) Sco X-l 
at a = 16 h 19 m 55.085 s , 5 = -15°38'24.90"; (2) 
NEIBC at a = 16 h 19 m 57.439 s , 6 = -15°37'24.0" 
(epoch J2000). 

2.2. Calibration 

The observations of 1504-166 were used to de- 
termine the delay offset and delay rate for each 
antenna, predominantly caused by the differences 
in the independent clocks and oscillators, and 
to check the general quality of the data. The 
delay and delay rates were determined by us- 
ing routine FRING in the AIPS software pack- 
age (Greisen 1988). The total flux density of 
the calibrator was measured to 2% accuracy us- 
ing the VLA observations of 3C286 to determine 
the absolute flux density scale of the observations. 
The milliarcsecond structure of 1504-166 was de- 
termined from the VLBA observations, and this 
model was used to determine the gain calibration 
for all of the data, with an accuracy of about 3%. 

The data for the VLBA and EVN observa- 
tions consisted of 8 independent streams, at 4 con- 
tiguous frequencies each with 2 polarizations, and 
the APT observations contained 2 data streams. 
After the 1504-166 calibrations, all data streams 
were combined to increase the signal-to-noise ra- 
tio (SNR) for the next calibration step using the 
relatively weak NEIBC. With the AIPS routine 
CALIB, we determined the residual phase error 
for each antenna every few minutes from the ob- 
servations correlated at the NEIBC position. At 
1.7 GHz, its correlated flux density was > 10 
mJy on the longest VLBA baseline; hence, an ac- 
curate antenna phase error could be determined 
from two minutes of data. This time is shorter 



3 



than the time-scale of phase changes except dur- 
ing rare periods of ionospheric turbulence. When 
the phase could not be connected unambiguously 
between consecutive phase solutions, the data be- 
tween these two times were omitted from the anal- 
ysis. 

At 5.0 GHz the calibration strategy was more 
complicated. The NEIBC correlated flux density 
decreased from 8 mJy at the shorter baselines to 3 
mjy for the longer baselines. At least two minutes 
(five minutes for the APT baselines) integration 
was needed to determine phase errors for these 
baselines. For low elevation observations or dur- 
ing periods of poor phase stability, phase solutions 
were ambiguous. Hence, at 5.0 GHz the NEIBC 
could reliably calibrate baselines shorter than 3000 
km. After this initial calibration of the shorter 
baselines, further self-calibration of Sco X-l suc- 
cessfully reinstated the longer baselines when Sco 
X-l was stronger than about 5 mJy. 

With the simultaneous observations of the 
small-diameter NEIBC, the calibration of Sco X-l 
allowed us to register all images over the four-year 
period on the same astrometric grid with an es- 
timated accuracy of < 0.1 mas (Bradshaw et al. 
1999). The image registration of Sco X-l was not 
altered when the data were further self-calibrated 
in order to add the longer spacings. This registra- 
tion assumed that the NEIBC was stationary with 
the same centroid at the two frequencies. This as- 
sumption is reasonable because: (1) The source 
had not varied at either frequency over the last 15 
years of VLA monitoring, (2) its structure is sym- 
metric, with an angular size of 2.5 mas at 5.0 GHz 
and 3.5 mas at 1.7 GHz, and (3) its spectral index 
is —0.4 (S oc v a ), suggesting the lack of a sig- 
nificant small-diameter opaque component, which 
might be variable. The source was not identified 
to a magnitude limit of 25 m (Malin, priv comm, 
cited in (Geldzahler et al. 1981)). 

2.3. Imaging, Self-Consistency and Pa- 
rameterization 

Deconvolved images of Sco X-l were made using 
the routine IMAGR in AIPS and with the Caltech 
software package Difmap (Shepherd 1997). Im- 
ages of the NEIBC were also made as a check 
on the quality of the calibration. Due to signif- 
icant variability of Sco X-l during the observa- 
tions, image artifacts associated with the aperture 



synthesis of variable sources often affected an im- 
age made from more than one hour of data. For 
the VLBA observations, when the SNR was suf- 
ficient (the flux density was more than about 5 
mJy), images could be made every hour at 1.7 
GHz from one 9-minute observation and at 5.0 
GHz from three 9-minute observations. If Sco X- 
1 was weaker than 5 mJy, longer periods of data 
were combined in order to obtain images with a 
useful SNR. 

The u-v coverage in the hour-to-hour VLBA 
snapshots was sparse, especially during the first 
and last hour of an observation, when the cov- 
erage contained large gaps. An rms noise of 
0.2 mJy/beam for each snapshot, a peak inten- 
sity in the range 1.5 and 20 mJy/beam, and a 
somewhat noisy phase calibration using the rela- 
tively weak NEIBC, when combined with the rel- 
atively poor u-v coverage and the source variabil- 
ity, conspired to make the determination of reli- 
able radio images ambiguous. For some snapshots, 
an unrestrained deconvolution of the images pro- 
duced physically unacceptable image characteris- 
tics which often had properties of the dirty beam 
which varied through the observation day. How- 
ever, the relatively simple structure of Sco X-l 
(a few bright small-diameter components along a 
well-defined position angle) made the deconvolu- 
tion less ambiguous than with a more complicated 
source. The justification of this assumption is that 
the use of this simple model (which means that 
only a small region of the image was searched for 
emission using the deconvolution algorithms) led 
to consistent results for the entire body of data. 

Two methods were used to judge the quality of 
the deconvolution of the images of Sco X-l. First, 
reasonable continuity of the images between con- 
secutive epochs was expected. Sudden changes in 
the source structure or image features which re- 
sembled the 'dirty-beam' suggested a poor de- 
convolution or bad data quality. The data were 
then scrutinized, the NEIBC calibration stabil- 
ity was checked, the source reimaged, and de- 
convolved. Secondly, the hourly data sets were 
also processed by the Difmap processing package. 
By fitting simple Gaussian component models di- 
rectly to the visibility data, some of the ambigu- 
ities of image-plane deconvolution were avoided. 
If the Difmap model did not agree well with the 
image, further checking was necessary. Finally, 
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the independent images made at 5.0 GHz and 1.7 
GHz had to be consistent; that is, the compo- 
nents at the two frequencies were nearly coinci- 
dent, although the 1.7 GHz lower resolution ob- 
servations did occasionally show more extended 
emission, as expected on physical grounds. The 
difference in angular scale between the 1.7 and 5.0 
GHz beam properties by a factor three also aided 
in determining real structure from beam-induced 
artifacts. Simulations were also made in order to 
quantify the reliability and proscribe the proper 
reduction methods. Data errors associated only 
with the expected level of stochastic receiver noise 
should cause little ambiguity in the images. Other 
potential data errors (unknown data drop-outs, 
strong source variations, rapid atmospheric phase 
variations over minute time-scales) were modeled, 
showing that relatively large errors could lead to 
significant image errors. For this reason we re- 
moved marginal data, especially for the longer 
baselines where these non-random errors occurred, 
with the concomitant loss of resolution for some 
snap-shots. 

For the 5.0 GHz EVN and APT observations, a 
similar imaging and reduction strategy was used. 
The NEIBC was only useful in calibrating the Eu- 
ropean baselines in the EVN and the Australian 
baselines for the APT. If Sco X-l was sufficiently 
strong (> 5 mJy for the EVN, > 10 mJy for 
the APT), then the other antennas could be self- 
calibrated using Sco X-l. However, the u-v cover- 
age for these arrays was significantly worse than 
that for the VLBA, so the data had to be averaged 
over several hours before reliable snapshot images 
could be made. For the APT observations, two im- 
ages separated by about three hours were made. 
For the EVN observation, three images, separated 
by about two hours, were made. 

The resolution of an image depended on the 
SNR and the elevation of the source. For VLBA 
observations at 1.7 GHz, all images were convolved 
to a full-width half-power resolution of 10 x 5 
mas at position angle 0°, although the resolution 
of the observations varied somewhat during the 
day. At 5.0 GHz, imaging using only the shorter 
spacings, produced the same resolution as that at 
1.7 GHz. These image pairs were used to deter- 
mine the spectral index distribution across Sco X- 
1. At 5.0 GHz, when the source was frequently 
stronger than 5 mJy and self-calibration of the 



longer baselines was possible, the resolution was 
about 4.5 x 1.5 mas during the middle part of the 
day. The resolution degraded to 6 x 2 mas at low 
elevation angles. For the occasional periods when 
Sco X-l was stronger than about 15 mJy, a reso- 
lution of 3 x 1 mas was obtained. 

The flux density, position and angular size for 
the discrete radio components in Sco X-l were 
determined in two ways: (1) fitting Gaussian- 
components directly to the image, and (2) fitting 
Gaussian-components to the visibility data. Al- 
though errors were estimated from each of the two 
methods, the difference in the radio source param- 
eters from the two analysis were used as a better 
accurate estimate of the real uncertainties. 

3. The Radio Images of Sco X-l 

3.1. 1999 June 11-13; MJD 51340-51342 

3.1.1. The Images and Component Parameters 

Figure 1 shows the changing structure of Sco 
X-l with 10 x 5 mas resolution at 1.7 GHz and 
5.0 GHz from the VLBA observations. These 
images, although only covering 40% of the ob- 
servation, form a consistent set from which gen- 
eral evolution properties are evident. The ver- 
tical dashed line shows the position of the bi- 
nary system, as determined from its radio parallax 
(Bradshaw et al. 1999). The skew dashed-line to 
the left shows the location of the component north- 
east of the core as a function of time. On MJD 
51342, the vertical line to the right shows the ap- 
proximate position of the component southwest of 
the core. The locations of the dashed lines on the 
1.7 GHz and 5.0 GHz columns are identical. 

More details of the changing structure of Sco 
X-l arc shown in Figure 2. These 5.0- GHz im- 
ages are at the highest resolution obtained during 
the 56-hour experiment. The image snapshots are 
separated by about 50 minutes for the VLBA ob- 
servations and about 2.5 hours for the APT and 
EVN observations. The APT and EVN images 
have a resolution of 10 x 5 mas, and most of the 
VLBA images have a resolution of 4.5 x 1.5 in po- 
sition angle 0°. The image field of view is 35 x 25 
mas for MJD 51340 and 51341 and 90 x 25 mas 
on MJD 51342. The vertical lines follow the same 
tracks as in Figure 1. 

Two estimates of the total flux density for Sco 
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X-l are shown in the top two plots of Figure 3. 
The flux density, shown in Figure 3a, was deter- 
mined from the VLA observations with determi- 
nation of the flux density made every 10 minutes. 
With a resolution of 3" at 5.0 GHz and 10" at 
1.7 GHz, these measurements contain the entire 
emission from Sco X-l. Figure 3b shows the to- 
tal flux density obtained from the VLB A images. 
The good agreement between the VLA and the 
VLBA measurements indicates that no significant 
large-scale structure associated with Sco X-l was 
resolved out by the VLBA observations. Figures 
3c, 3d and 3e show the flux density for each of the 
three components. In Figure 4, the spectral in- 
dex a of each of the three components during the 
dual-frequency VLBA observations is shown. 

3.1.2. A Description of the Changes in Sco X-l 

The 1999 June observations form a basis for the 
interpretation of the radio emission from Sco X-l. 
Thus, we will describe the evolution seen in the 
images in some detail. 

The basic morphology of Sco X-l is simple. The 
source is usually composed of three relatively com- 
pact components: a radio core that is nearly co- 
incident with the binary system, a compact NE 
component moving away from the core; a weak 
SW component moving away on the opposite side 
of the core. The three components lie along an 
axis of position angle 54° to an accuracy of a few 
degrees. 

Changes associated with the emission from Sco 
X-l are as follows: 

• MJD 51340.1 to 51340.4. The NE com- 
ponent reached a maximum intensity, with 
a = -0.5, at MJD 51340.2. Subsequently, it 
faded, became more extended with a steep- 
ening spectral index, and fell below the 
VLBA detection level at MJD 51340.4. It 
did not reappear over the next ten hours in 
the images from the APT and EVN. Dur- 
ing this period, the radio core flux density 
steadily increased with a > 0.0. 

• MJD 51340.4 to 51341.3. The radio core 
decreased in flux density from 20 mJy to 8 
mjy in a few hours (the rms error for a flux 
density measurement is about 0.5 mJy.) No 
emission to the north-east was present. At 



MJD 51340.7 the core flux density began to 
increase, and by MJD 51341.0, the flux den- 
sity had reached 20 mJy. The images clearly 
show that this rise in flux density was asso- 
ciated with a component emerging from the 
core to the northeast. Between MJD 51341.0 
to 51341.3, this new NE lobe moved outward 
and decreased in flux density. Emission be- 
tween the lobe and core was also present. 
The core flux density generally decreased in 
flux density during this period, but a small 
flare occurred at MJD 51341.3. 

• MJD 51341.3 to 51341.8. After moving 
about 20 mas from the core, the NE com- 
ponent brightened, slowly at first, and then 
abruptly flared from 2.8 mJy to 20 mJy in 
less than one hour. The spectrum quickly 
flattened to a = —0.5, then NE flux den- 
sity slowly decayed to 1.5 mJy over the next 
15 hours and its spectral index dropping to 
a w —1.2. The radio core, during this pe- 
riod, remained relatively constant at 4 mJy, 
except for a flare to 8 mJy at MJD 51342.3 

• MJD 51341.8 to 51342.2. The radio core 
brightened from 4 mJy to 20 mJy, reach- 
ing the peak at MJD 51342.0. The radio 
core then decreased in flux density and by 
MJD 51342.1 extended emission towards the 
northeast was clearly present. However, un- 
like the core flare and subsequent expulsion 
of emission to the NE about 1.3 days earlier, 
this extended emission decreased below the 
detection level by MJD 51342.2 and with no 
apparent motion. The SW component, seen 
better in the lower resolution 1.7 GHz im- 
ages, is visible during this period. The NE 
component decreased slowly in flux density. 

• MJD 51342.2 to 51342.4. For the last 
four hours of the observation, the NE com- 
ponent brightened from 1.5 to 4 mJy and the 
spectral index flattened again to —0.5. This 
NE component flare is less spectacular than 
the one about 0.8 days earlier, but has sim- 
ilar properties. The core remained steady 
at about 4 mJy with no resolvable extended 
structure. 
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3.2. 1998 February 27-28; MJD 50871-2 

The only previous Sco X-l observation on con- 
secutive days occurred on 1998 February 27-28, 
from two six-hour VLBA observations separated 
by 18 hours. The upper part of Figure 5 shows 
three sets of snapshots, each 52 minutes apart 
on MJD 50871. The lower two panes show the 
structure on the following day when the source 
was weaker. Figure 6 shows the flux density of 
the three components on both days. Although 
there was some indication of extended structure 
between the core and the NE component, the high- 
resolution 5.0 GHz images show that most of the 
emission was contained in compact components. 

On February 27, the core spectrum was inverted 
with flux density about 10 mJy at 5.0 GHz. The 
NE component had a flux density > 15 mJy and 
was moving from the core in the same position 
angle as that in the 1999 June observations. Near 
the end of the observation period the flux density 
decreased, its spectral index became steeper and 
the angular size increased. This behavior will be 
discussed in more detail in §5.7 and with Figure 
14. 

Sco X-l was much fainter on the following day 
and individual snapshots were too noisy to show 
rapid changes. At the bottom of Figure 5, the 1.7 
GHz and 5.0 GHz images, both with 10x5 mas res- 
olution, were made from the entire 5-hour obser- 
vation. The core was unresolved with a peak flux 
density of 3.0 mJy, and decreased during the day. 
The SW component, clearly detected at about 2 
mJy, was also relatively compact. The radio emis- 
sion northeast of the core was complex but con- 
fined along the major axis of the radio structure, 
between 30 mas to 60 mas from the core at 1.7 
GHz. In the 5.0 GHz image, only two features in 
this northeast emission region were detected above 
the noise level. A linear extrapolation of the NE 
component motion on MJD 50871 intercepts the 
most distant component of emission in the north- 
east direction on MJD 50872. Motion of the SW 
component was detected on both two days, and is 
discussed in the following section. 

3.3. The Component Velocities and Orien- 
tations 

The 1998 February and 1999 June images show 
the major dynamical properties of Sco X-l. A 



compact component was always present near the 
binary system. It varied significantly in flux den- 
sity and occasionally faint emission extended sev- 
eral mas towards the north-cast. A compact com- 
ponent was usually detected northeast of the radio 
core and its motion was clearly directed away from 
the radio core. Its flux density also varied signif- 
icantly over hour time scales. A compact com- 
ponent southwest of the core was detected about 
50% of the time and was on average about 10 times 
weaker than the NE component. 

The relative simplicity of these images of Sco 
X-l is contrasted by the more complicated emis- 
sion from other well-studied Galactic-jet radio 
sources (Mirabel & Rodriguez 1999). Although 
their structures tend to lie on a well-defined axis, 
they are often composed of many components on 
both sides of the radio core. For Sco X-l, we find 
that only one compact component was detected on 
the NE and SW sides of the radio core (See Figure 
11 for a possible exception). 

The changing separation of the NE and SW 
components from the core for the 1999 June (top) 
and 1998 February (bottom) observations is shown 
in Figure 7. For the 1999 June observations, the 
velocity in the plane of the sky for the NE com- 
ponent #1, detected for 6 hours on MJD 51340, 
was v = 0.73 ± 0.07 mas hr" 1 = 0.28c 1 . The 
NE component #2 was not resolved from the 
core until about MJD 51341.0. Then, from MJD 
51341.0 to 51341.35, the NE component speed was 
v = 1.74 ± 0.16 mas hr" 1 =0.68c. After MJD 
51341.4, when the NE component flared, the com- 
ponent speed decreased to v = 1.25 ± 0.05 mas 
hr _1 = 0.49c. Although the NE component var- 
ied considerably in flux density during the last 
24 hours of the observations, its speed remained 
nearly constant. 

For the 1998 February observations, the ve- 
locity of the NE component on MJD 50871 was 
v = 1.11 ± 0.06 mas hr _1 = 0.43c over the five 
hours of observations, with no significant depar- 
ture from linearity except for the last frame. The 
single point for the NE component on MJD 50872 
was taken as the position of the most distance part 
of the NE component seen in Figure 5. 

Motion of the SW components was more dif- 



For a measured distance to Sco X-l of 2.8 kpc, one mas = 
2.8 AU = 4.19 X 10 13 cm. A speed of 1 mas hr" 1 = 0.388c. 
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ficult to detect. On MJD 50872 the component 
was both strong and long lasting and moving 
away from the core with v = 0.6 ± 0.2 mas hr _1 . 
On MJD 51342 the velocity was v = 0.5 ± 0.3 
mas hr _1 . In Figure 7, the dashed-lines passing 
through the positions of the SW components are 
not linear fits to its position, but to 50% of the 
fitted angular separation to the NE compo- 
nent (Figure 12 shows a plot of this ratio from 
all Sco X-l observations). Thus, we surmise that 
the SW component was also moving away from 
the core, but at a projected velocity of about 50% 
of that of the NE component. These motions are 
described in more detail in §4.1. 

The orientation of the NE and SW components 
with respect to the core for the 1999 June and 1998 
February observations are shown in Figures 8a and 
8b. The radio structure of Sco X-l remained near 
position angle 54° during these these two obser- 
vations, and the component orientations from the 
previous observations of Sco X-l also lie within a 
few degrees of 54°. Hence, any long-term variation 
in the orientation of the radio emission from Sco 
X-l is less than about three degrees. It is possible 
to fit the position angle distribution in 1999 June 
observations with a mean value of 54° plus a sinu- 
soidal term of « 3° with a period of between two 
to four days. Modeling of the binary system X- 
ray QPO's suggest that the neutron star/accretion 
disk interaction could produce a 5° precession over 
several days (Titarchuk & Oshcrovich 2000). 

3.4. The 1995-1996 VLB A Observations at 
5.0 GHz 

The first four VLBA observations of Sco X- 
1 were made on 1995 August 19 (MJD 49948), 
1996 March 16 (MJD 50158), 1996 September 16 
(MJD 50340) and 1996 August 3 (MJD 50663), 
and contour plots of the radio emission are shown 
in Figure 9. These observations had lower sensi- 
tivity than later observations, were only at 5 GHz, 
and were about four hours in duration. The radio 
emission from Sco X-l during these observations 
was relatively weak so that images made from one- 
hour slices of the data were too noisy in order to 
determine variability and motion. Nevertheless, 
all four observations showed similar properties in 
the emission of Sco X-l. 



• MJD 49948. The radio core had a peak 
flux density of 0.32 mJy. There was a trace 
of both the NE component and the SW com- 
ponent 19 mas and 12.2 mas from the core 
respectively, each with a peak flux density of 
< 0.20 mJy. 

• MJD 50158. The radio core had a peak 
flux density of 0.50 mJy. There was a slight 
indication of the NE component 14 mas from 
the core, with a peak flux density < 0.25 
mJy. 

• MJD 50340. The source was stronger dur- 
ing this observation and all three compo- 
nents were present well above the noise. The 
SW component had an integrated flux den- 
sity of 2.5 mJy, located 5 mas from the core. 
The core flux density increased from 0.30 
mJy to 0.60 mJy during the four-hour ob- 
servations. The NE component was about 8 
mas from the core with a peak flux density 
of about 0.6 mJy. The distance ratio of the 
NE and SW components from the core was 
about 2:1. 

• MJD 50663. The core had a peak flux den- 
sity of 2.0 mJy and the NE component was 
detected, about 26 mas away, with a peak 
flux density of 0.35 mJy. Any SW compo- 
nent was < 0.15 mJy. 

The relative weakness of Sco X-l in 1995- 
1996, compared with later observations, is con- 
sistent with the general characteristics of the 
long-term flux density behavior of Sco X-l, 
as shown from the Green Bank Interferome- 
ter monitoring between MJD 50600 and 51700 
(see www.gb .nrao.edu/ fgdocs / gbi /pubgbi/ ScoX- 
1). At the monitoring frequency of 8 GHz, the 
quiescent level of Sco X-l remained at 5 mJy for a 
period of 100 to 150 days, but there were periods 
of enhanced emission at 10 to 15 mJy which per- 
sisted for 50 to 100 days. About 5% of the time 
there were large outbursts (> 50 mJy), which 
lasted over a few days. These outburst occurred 
randomly and were not associated with the peri- 
ods of enhanced emission from Sco X-l. 

3.5. 1998 August 29-30, MJD 51054-5 

Figure 10 displays the structure of Sco X-l for a 
five-hour observation on August 29-30, 1998. Be- 
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cause of the weakness of the source, each image 
contains data from the entire five-hour observa- 
tion. All three components were present, although 
the core was barely detectable at either frequency. 
The NE component was clearly extended and con- 
tained several emission peaks along the source ma- 
jor axis. The SW component was compact and 
relatively strong. The appearance of Sco X-l on 
this day was similar to that on MJD 50872, shown 
at the bottom of Figure 5. However, the linear 
size of the source on MJD 51054 was a factor of 
three smaller. Again, the ratio of the separation 
from the core of the NE and SW components was 
about 2:1. 

3.6. 1997 August 22, MJD 50682 

During this five hour VLBA observation, Sco 
X-l was extremely bright and nearly all of the 
emission was contained within 5 mas of the core, 
as shown in Figure 11. The high-resolution image 
showed emission moving outward from the binary 
system, particularly on the north-east side of the 
core. The relatively bright spot of emission at the 
extreme edge of the extended emission was ad- 
vancing at a velocity of 1.4 ± 0.3 mas h _1 . Some 
emission was also detected SW of the core. This 
morphology and evolution was similar to that seen 
on MJD 51341.0 when the second NE component 
was ejected from the core. The spectral index of 
the emission within 5 mas of the core rose from 
a = to a — 0.3 during the five hour observation. 

A faint component, about 20 mas NE of the 
core, was also detected. This is the only example 
of two radio components on one side of the core. 
However, the distant component is extremely weak 
with no measured motion. The radio structure on 
this day will be discussed in more detail elsewhere 
(Geldzahler et al. 2002). 

4. Orientation of the Binary system 

4.1. Doppler-Beaming of NE and SW 
Components 

The kinematic results discussed in the section, 
have already been reported in Paper I. We present 
more detail in this section. 

The images from the 1998 February and 1999 
June observations showed that there was often a 
compact NE component moving away from the 



core. Although this component varied signifi- 
cantly in flux density, its velocity remained rela- 
tively constant over periods of many hours. When 
the NE component was well-separated from the 
core, its measured speeds in the plane of the sky 
were 0.28c, 0.43c, 0.49c and 0.68c. On MJD 51341 
and MJD 50682 a NE component emerged from 
the core. Because of a lack of sufficient resolu- 
tion, the emergent velocity was not well-defined, 
although it was about 0.5c. 

Motion of the SW component away from the 
core was observed during the 1998 February and 
1999 June observations, with v = 0.6 ± 0.2 mas 
hr _1 , and v = 0.5 ± 0.3 mas hr _1 , respectively. 
Further detection of the motion was not possible 
because of the weakness and intermittent detec- 
tion of this component. A comparison of the NE- 
to-corc and the SW-to-core component distance 
ratio of all the observations is shown in Figure 

12. There is no systematic variation with distance 
and the average ratio is 0.51 ± 0.02. Since the NE 
component was moving away from the core, we in- 
fer that the SW component was also moving away 
from the core in the opposite direction with 51% 
of the velocity of the NE component. 

Another comparison of the radiative properties 
of the NE and SW components is given in Figure 

13, which displays the spectral index versus flux 
density for the two components in Sco X-l from all 
of the observations whenever it was determined. 
For comparison, the spectral properties for the 
core are shown (Geldzahler et al. 2002). When 
the flux density of the NE component was stronger 
than 4 mJy its spectral index was a ~ —0.6 but 
decreases to < —0.9 at fainter flux densities. The 
SW component was detected less than 50% of the 
time and often only at 1.7 GHz. It greatest flux 
density was 5 mJy at 1.7 GHz and the spectral 
index was about —0.6; hence, the flux density at 
5.0 GHz is about half of that at 1.7 GHz. We 
conclude that on average the apparent brightness 
of the SW component is about ten times fainter 
than that of the NE component. Furthermore, the 
spectral index of the SW component at its bright- 
est was about equal to that of the NE component 
at its brightest. We will discuss the flux density 
ratio between the two lobes in more detail in §5.9 

The relative positions and velocities of the NE 
and SW components, their relative flux density (in 
a statistical sense since both components vary con- 
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sidcrably), and their spectral properties suggest 
strongly that the NE and SW components have 
similar properties, with the observed emission dif- 
ferences resulting from relativistic aberration and 
Doppler beaming. The observed apparent motion 
in the plane of the sky of a component moving 
with a velocity v (usually given in term of c as 
=v/c), at an angle of 9 to the line of sight is (1): 



0ne — 



0sw 



sin(6>) 
1 - cos(0) 
sin(0) 



(1) 



(2) 



1 + cos(6») 

where 0ne is the observed velocity in the plane 
of the sky of the component approaching the ob- 
server and 0sw is the velocity of the receding com- 
ponents. 

We have measured a range of speeds for the 
NE component. However, the nearly constant po- 
sition angle of the axis of Sco X-l in the plane 
of the sky (Figure 8) strongly suggests that the 
direction of the space velocity of the lobes with 
respect to the line of sight is also constant. To 
determine an average speed of the two Sco X-l 
components, we considered the time ranges MJD 
51341.0 to 51342.4 and MJD 50871.3 to 50872.4. 
These periods are instances when the NE compo- 
nent speed was determined over a long period of 
time and the SW component was relatively strong. 
These two periods yield a measured proper motion 
oiv = 1.25 ± 0.05 and v = 1.11 ± 0.06 mas hr" 1 . 
We will thus assume that the average speed of the 
NE component of Sco X-l of 1.18 ±0.08 mas hr" 1 . 
At a distance of 2.8 ±0.3 kpc for Sco X-l, the aver- 
age speed of the NE component in the plane of the 
sky is then 0ne = 0.46 ± 0.08, with the distance 
uncertainty incorporated in the error. Next, we 
adopt the ratio of the SW speed to the NE speed 
as 0.51±0.02. This ratio is a function of the speed, 
but should be the appropriate value that we have 
adopted. Equations (1) and (2) become: 



0sw 



= 0.51 ±0.02 



1 - cos(0) 
1 + cos(0) ' 

sin(0) 
1 - cos(0) ' 

= 0.45 ±0.03; = 44° ±6°. 



0NE 

NE = 0.46 ±0.08 = 
from which we obtain 



(3) 
(4) 

(5) 



The quoted errors are the one-sigma uncertainties 
and include the measurement errors. 

The determination of the Doppler boosting as- 
sociated with the flux densities of the NE and SW 
components requires a specific model of their inter- 
nal structure, particle motion and magnetic field 
geometry as discussed in §5.5. If we assume that 
the moving components are transparent, with all 
particles moving at the bulk velocity and radiating 
isotropically, then the flux density ratio, R is given 
by (Blandford & Konigl 1979; Cawthorne 1991): 



R 



-i k — a 



SW 



0. 



NE 



(6) 



where k is a geometric factor which equals 2 for 
an optically thin source, and a is the spectral in- 
dex of the source emission. From the measured 
a = —0.6 we obtain a predicted flux density ratio 
of R=0.17, which is somewhat larger than that ob- 
served (see §5.9 for more details). Thus, the space 
motion derived solely from the relative motions of 
the two components is in reasonable agreement, 
considering many uncertainties about the radia- 
tion properties and variability of the lobes. 

4.2. Sco X-l Binary Masses 

Assuming that the radio beam lies along the 
rotation axis of the accretion disk which in turn 
lies in or close to the orbital plane of the binary 
system, we can better determine some of the prop- 
erties of the binary. The relative constancy of the 
radio axis over five years are supports the above 
association. 

The mass function for the Sco X-l binary 
has been derived from the measured radial ve- 
locity of the Hell A-4686 line, presumably from 
the accretion disk, and is given by F(M) = 
(M 2 sin i) 3 /(M 1 + M 2 ) 2 = 0.016 ± 0.004 (our 
estimated error), where Mi and M 2 arc the 
masses of the neutron and secondary star, re- 
spectively, and i is the inclination of the orbit 
(Cowley & Crampton 1975). The binary is com- 
posed of a degenerate star, almost certainly a neu- 
tron star (Hasinger & van der Klis 1989). Mea- 
surements of many neutron star systems pro- 
duced estimated masses between 1.2 to 1.6 Mq 
(Lewin et al. 1995). The intrinsic X-ray luminos- 
ity of Sco X-l, for a distance of 2.8 kpc, also 
suggests a neutron star mass of about 1.4 Mq. 
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The calculated secondary mass M2, from the best 
value of the mass function, inclination and neutron 
star mass, is (0.63 ± 0.26) Mq. Such a low com- 
panion mass is a problem for most models of the 
energetics in Sco X-l. First, the star may not fill 
its Roche lobe resulting in an accretion rate that is 
too low; perhaps the companion is an evolved star. 
Secondly, infra-red spectroscopic observations of 
the binary system suggest that the secondary star 
is earlier than G5 (Bandyopadhyay et al. 1999). 
However, the X-rays from Sco X-l in all likeli- 
hood heat up the atmosphere of the secondary, 
puffing out its atmosphere until it overflows the 
Roche lobe reulsting in mass accretion onto the 
compact star. 

The secondary star mass estimates can be in- 
creased in the several ways. If the degenerate star 
mass is > 3.5 Mq, then the secondary star would 
have a mass > 1 M . But, then what is the nature 
of the degenerate star? A mass function > 0.025 
would also increase the secondary star mass associ- 
ated with a 1.4 M Q neutron star. This implies that 
the optical lines are coming from a much larger 
region than the accretion disk near the neutron 
star. Thus, the masses of the neutron star and 
the secondary member of the Sco X-l binary are 
still uncertain even with knowledge of the orbit 
inclination. 

5. The Nature of the NE Component 
5.1. Angular and Linear Size 

The angular size of the NE component could 
be measured with an accuracy of 0.5 mas when it 
had a flux density > 5 mJy. Figure 14 shows a plot 
of the angular size and motion of the NE and core 
components for the seven snapshots on MJD 50871 
(see Figure 5 also) when both the core and NE 
components were relatively strong. The angular 
sizes were determined by fitting the visibility data 
directly to an elliptical Gaussian model and by 
fitting the component brightness on the image and 
deconvolving the effects of the beam shape. Both 
methods gave consistent results. 

The comparison of the positions and sizes of the 
core and the NE components provides a good in- 
dication of the robustness of these measurements. 
The core was about 10 mJy, relatively stationary 
and extended in the direction along the source 
major axis, as was commonly observed at other 



epochs. Alternatively, the NE component was ex- 
tended orthogonal to the direction of motion and 
was clearly moving away from the core. Since the 
effects of noise and poor phase calibration will blur 
all of the emission in the same way, the different 
size and orientation of the core and NE compo- 
nents strongly suggests that these measured an- 
gular sizes and orientations are valid. 

For snapshots 1 through 3 in Figure 14, the 
measured average angular size and orientation of 
the NE component itself was (1.5 ± 0.15) x (0.9 ± 
0.2) mas in position angle (pa) 140° ± 15°. The 
angular size in snapshot 4 is marginally bigger and 
the size for snapshot 5 is 2.4 x 1.3 mas in pa 120°. 
In the 50-minute separation between snapshot 5 
and 6, the component increased in area by a fac- 
tor of 4.5, to 5.0 x 2.9 mas in pa 150°, an expan- 
sion at a velocity close to c. In snapshot 7 the 
size of the component did not increase substan- 
tially, but the bulk motion of the emission may 
have decreased. The angular size of the NE com- 
ponent was also accurately determined on MJD 
51341.4 and MJD 51342.3. The angular sizes were 
(1.6 ± 0.1) x (0.7 ± 0.3) mas in pa 135° ± 10 and 
(2.4 ± 0.3) x (0.8 ± 0.4) ma in pa 155° ± 30°, re- 
spectively, in good agreement with the minimum 
size of the NE component on MJD 50871. 

The NE component is clearly extended in the 
direction perpendicular to its motion. In addi- 
tion, the measured size in this direction of mo- 
tion is slightly increased by the component mo- 
tion during each snapshot over a 33-minute span. 
During this period, the component moves 0.6 mas; 
hence, the measured size of 0.9 mas corresponds 
to a true size of about 0.7 mas after removing the 
blurring from the motion. Thus, wc derive a mini- 
mum component size (full-width at half power) of 
1.5 x 0.7 mas, corresponding to 6.3 x 10 8 km by 
2.9 x 10 8 km with the major axis perpendicular to 
the orientation of the source. More detailed mea- 
surement of the NE component structure, other 
than its angular size, cannot be obtained with the 
present resolution and SNR. 

The flux densities and spectral indices for the 
NE components are also listed in the Figure 14. 
The relationship between the angular size, flux 
density and spectrum of the NE component are 
discussed in §5.7. 
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5.2. Component Energetics 

The calculation of the energetics of the NE com- 
ponent used the following parameters: a source 
distance of 2.8 kpc (8.6 x 10 16 km), a flux den- 
sity at 1.7 GHz of 30 mJy, a = —0.6, with the 
radio cutoff frequencies of 10 7 to 10 11 Hz. The 
volume of the component (ellipsoid of diameters 
1.5x1.5x0.7 mas with the narrow dimension in the 
direction of the component advance) is 6.1 x 10 40 
cm 3 . The integrated flux density is 6.5 x 10~ 21 
erg s^cm -2 , the total observed radio luminosity 
is 6.1 x 10 30 erg s _1 . We assume a factor of 3 to 
correct for the Doppler boosting of the NE com- 
ponent (and Doppler attenuation of 3 for the SW 
component) as implied by the average flux density 
ratio of the two component . The radio luminosity 
in the frame of the source is 2.0 x 10 30 erg s -1 , cor- 
responding to a brightness temperature of about 
1 x 10 9 K. The total energy in relativistic parti- 
cles E e , and total energy in the magnetic field Ejj, 
assuming the above radiating volume, are 

£ e (ergs) = 4.0 x 10 37 H~ 15 (7) 

E H (ergs) = 2.4 x 10 39 H 2 , (8) 

where H is the magnetic field in Gauss. 

The radiative lifetime of the emitting particles 
is about one year if magnetic field energy and ki- 
netic energies are about equal. For example, if we 
assume that the protons contain 100 times the en- 
ergy as the electrons, the equipartition magnetic 
field is about 1.0 G with a minimum total energy 
of 6 x 10 39 erg in the component. If the electron 
energies dominate, then the minimum total energy 
is about 5 x 10 38 erg and has a field strength of 
about 0.3 G. Such a long radiative lifetime is not 
possible for the NE component since it is variable 
on an hour time-scale and a size of 5 x 10 8 km 
(30-min light travel time). This problem will be 
discussed in §5.6. 

5.3. What is the NE component? 

Emission confined to relatively discrete com- 
ponents, and which show relativistic motion have 
been observed in many Galactic X-ray binary sys- 
tems (Mirabel & Rodriguez 1999). Three expla- 
nations are generally given for the energetics of 
these components: (1) Ejected clouds, where the 
components are clouds of radiating particles, ex- 
pelled from the binary in a preferred direction, 



but containing their own energy source; (2) Inter- 
action within beams; where the energy flow within 
twin-beams is collimated by an accretion disk of 
the massive object and interacts with embedded 
matter or entrained material to produce radiat- 
ing electrons; (3) a working surface where the en- 
ergy flow in the efficient twin-beams is relatively 
invisible until the flow impinges on external ma- 
terial and forms a bright, small region of highly 
relativistic particles. The beam also penetrates 
through the external medium at relativistic veloc- 
ities. We suggest that the moving components on 
opposite sides of Sco X-l are regions of intense ra- 
dio emission, generally denoted as lobes, formed 
from the interaction of a twin-beam with the am- 
bient medium in a confined region called a working 
surface. 

There are two properties of the components 
in Sco X-l that that are best explained by 
this lobe model. First, the rapid variability 
and lack of systematic decay of these compo- 
nents implies that the radiating electrons have 
a lifetime of less than one hour and are con- 
stantly resupplicd with energy. In fact, the NE 
component variability is similar to that of the 
core. In contrast, the emission from the rela- 
tivistic moving components associated with the 
other well-studied Galactic-jet radio sources are 
not as variable and tend to decrease fairly uni- 
formly as they move away from the radio core 
(Fender et al. 1999; Hjellming & Rupen 1995). 
Secondly, the observations of Sco X-l suggest that 
only one pair of components (NE and SW) occur 
at any given time (see Figure 11 for a possible 
exception). The core does show activity and ex- 
tended emission towards the northeast direction, 
but more than one compact component is not ob- 
served on either side of the core. For the other 
Galactic-jet sources, several components are usu- 
ally visible on one side of the core and they are of- 
ten associated with X-ray and radio flaring events 
in the core region. 

However, the later stages of the lobe evolution 
in Sco X-l did show a more complicated structure. 
For example, the radio emission from Sco X-l on 
MJD 50872 and MJD 51054 had an extended NE 
component containing several bright regions. On 
MJD 50872 this extended NE component is clearly 
associated with the moving compact component 
observed on the previous day and consists of an 
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extended lobe towards the core with a few small 
brighter regions. Perhaps this extended emission 
is back flow from the lobe region seen in many 
extragalactic objects (Black et al. 1992). 

5.4. The Advance Velocity of the Lobe 

The speed of advance for three different lobe- 
pairs were measured over the range 0.32< (3 <0.57. 
The NE lobe on MJD 51340 moved at 0.32c just 
before evaporating. The lobe on MJD 50871 
moved at 0.43c. The lobe on MJD 51341-2 moved 
initially at 0.57c, but abruptly decreased to 0.46c 
just after a lobe flare. The advance speed during 
each epoch was remarkably constant over periods 
of many hours, up to one day. 

A prediction of the advance speed of the lobe 
requires a detailed model of the interaction of the 
beam with the ambient medium, A more general 
approach (Blandford & Konigl 1979) has derived 
the velocity of the advance of a lobe under the as- 
sumption that the relativistic shock produces an 
ultra-relativistic equation of state on both sides of 
the shock. This region is generally called the work- 
ing surface and this region produces a radio hot- 
spot which may be a small, but intense, part of 
the lobe. If p\ and p 2 are the pressures and v\ 
and V2 are the velocities in the unshocked (beam) 
and shocked (working-surface) regions, then 

Vl_ _ Pl+ 3P2 V2_ _ / 3pi + p 2 , . 

c ^y3(3pi+p 2 ) c ~ y 3( Pl +3p 2 )' ( ' 

When the beam pressure is much greater than the 
shock pressure, the advance speed is 0.33c. For 
more dominant shock pressure, the advance speed 
eventually reaches c. This model gives results 
which are consistent with the range of observed 
speeds using reasonable pressure ratios across the 
hot-spot. However, the small variation of the ob- 
served speed with time is still puzzling. 

5.5. Lobe Emission Geometry 

The ultra-relativistic fluid produced at the 
working surface fills the larger volume of the ra- 
diating lobe. Figure 15 shows a schematic of the 
general appearance of the lobe based on some very 
simple assumptions. There are two velocities asso- 
ciated with the lobe: the average bulk motion ve- 
locity of « 0.45c, and a diffusion velocity of 0.57c 
associated with this relativistic gas. We assume 



further that the radiative lifetime of the emitting 
particles is about 30-minutes as required by the 
variability time-scale. The 'wagon-wheel' diame- 
ter (in the frame of the lobe) is then 6.2 x 10 8 km 
in diameter, and the thickness is the distance trav- 
eled by the lobe in 30 minutes, 2.4 x 10® km. The 
lobe, however, viewed at an angle of 44° to the 
line of sight, would make the component appear 
somewhat more circular than the above model. 

A deeper analysis of the physical processes in 
the lobes requires detailed hydrodynamic jet mod- 
eling and the radiative transfer of the electrons. 
The modeling of (Mioduszewski et al. 1997), for 
example, shows striking similarities between the 
structure of Sco X-l with some of their simulated 
models at a viewing angle of 45° to the line of 
sight (see their Figures Id, 5c, 16a). The pre- 
dicted peak brightness ratio between receding and 
advancing lobes is R=0.05. However, this model 
assumes an isotropic magnetic field distribution. 
In any case, there seems to be a small family of 
models whose parameters can be chosen to match 
the morphology as well as the physical properties 
and characteristics of the lobes. 

5.6. Lobe Emission Lifetime 

The lifetime of the radiating plasma in the lobe 
cannot be much longer than 30 minutes, the vari- 
ability time-scale. The two most likely sources of 
energy depletion are from synchrotron losses or 
from adiabatic expansion of the plasma. First, if 
equipartition holds between the magnetic energy 
and particle energies as discussed in §5.2, then the 
time scale for synchrotron losses is about one year 
for the equipartition field strength of about 0.5 
G. The radiative lifetime decreases with as H~ 15 ; 
hence a one hour decay time requires a magnetic 
field strength of H=300 G. The derived lobe ener- 
gies are then 

E H (crgs) = 2 x 10 44 for H= 300 G (10) 

^e(ergs) = 7 x 10 33 for H= 300 G. (11) 

The magnetic energy is then a factor of 10 8 larger 
than the equipartition value, even assuming a den- 
sity of protons about 100 times that of the elec- 
trons. Magnetic field amplification just in advance 
of the shock has been suggested (De Young 1980; 
Hughes & Allen 1981), although not to this ex- 
treme. With a lobe brightness temperature < 
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10 9 K, the rate from inverse-Compton losses is 
much less than for synchrotron losses even with 
this large magnetic field. Such a large magnetic 
field requires a lobe energy density which is about 
10 5 times larger than the minimum energy value. 
The outer part of the radio lobe would also have 
a steeper radio spectrum than near the working 
surface and the 1.7 GHz angular size would be 
larger than that at 5.0 GHz. Unfortunately, the 
present observations just barely resolve the lobe 
at 5.0 GHz and such a spectral gradient is not 
observable. 

The radiating electrons generated within the 
working surface can also rapidly lose energy if 
they expand adiabatically to fill the lobe region. 
Assuming that the radio emission is optically 
thin (the radio spectral index a rts —0.6 is con- 
sistent with this) and neglecting the motion of 
the lobe for the moment, the steady-state emis- 
sion profile of an expanding spherical cloud which 
contains a steady source of accelerated electrons 
within a radius ro, the working surface region, is 
(Ball & Vlassis 1993; Hjcllming & Han 1995) 

S(r) = ^(r/ro) 4 "- 1 (12) 

where So is related to the inserted energy per unit 
time in the working surface volume, r is the dis- 
tance from the center of the working surface. For 
a = —0.6, the intensity at r/ro = 1.22, is already 
less than half of the intensity at the working sur- 
face. Hence the adiabatically expanding electrons 
need not diffuse far from the working surface in 
order to lose most of their energy. From the mea- 
sured lobe size of about 1 mas and a diffusion 
speed of 0.57c, the electrons will lose most of their 
energy in about 40 minutes. The bulk motion of 
the lobe is somewhat less than the diffusion speed 
and does not substantially change the above calcu- 
lation. Adiabatic expansion does not modify the 
spectral index of the emission across the lobe, in 
contrast to the spectral steeping with synchrotron 
ageing. In the case of pure adiabatic expansion, 
the working surface is an appreciable part of the 
observed radio lobe. 

A combination of both loss mechanisms may 
be operating in the lobes of Sco X-l. For exam- 
ple, the volume close to the working surface may 
have a large magnetic field which completely dom- 
inates the plasma dynamics. Further away from 
the working surface, adiabatic expansion may be 



the dominent loss mechanism. This would reduce 
the total energy content needed to supply the lobe 
emission compared with that needed if only syn- 
chrotron losses were operating. 

5.7. Lobe Expansion 

A slow expansion of the lobe angular size over 
an hour or two time-scale was observed several 
times during the observations of Sco X-l. The 
best example of the expansion of a lobe was on 
MJD 50871 after frame 3 in Figure 14. Between 
frames 5 and 6 the angular size increased by a 
linear factor of 2.2 in 50 minutes. As discussed 
in the above section, for pure adiabatic expansion 
with a uniform input energy source, the flux den- 
sity in frame 6 should have decreased (2.2) 4 " -1 
(Hjellming & Johnston 1988) or about 4% of that 
in frame 5. However, the flux density at snapshot 
6 was 60% of that in snapshot 5. 

Without higher resolution observations of the 
lobe we can only speculate on this slow lobe ex- 
pansion phenomenon. Since adiabatic expansion 
alone would produce a much more rapid decrease 
of flux density than that observed, we suggest that 
the increasing size and decreasing intensity of the 
lobe reflects the energy flow and size of the work- 
ing surface. The somewhat decreasing spectral in- 
dex could be related to the acceleration process in 
the working surface, or the effect of synchrotron 
ageing. Note that this lobe, although decreas- 
ing in flux density at the end of MJD 50871, is 
present on MJD 50872 as an extended lobe along 
the axis of the source. Hence, this lobe was rein- 
vigorated sometime after the end of the observa- 
tions on MJD 50871. Much more data are needed 
in order to explain the general appearance and 
evolution of these components. 

A similar expansion and spectral steepening of 
a lobe occurred during the first four hours on MJD 
51340. From the images in Figure 2, the compo- 
nent clearly expanded, the flux density decreased 
and the spectral index steepened to —1.5 until the 
component fell below the detection level of the ob- 
servations. 

5.8. Energy Flow Velocity in the Beam 

Paper I we showed that the the flux density 
variations of the radio components have good tem- 
poral correlation when assuming a model of energy 
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flow from the core to the lobe where: (1) a core 
flare is associated with an event near the binary 
system (2) a burst of energy associated with this 
event travels down the twin-beam at a velocity of 
(3j , (3) this increased energy flux intercepts the NE 
and SW lobes and produces a flare. 

The appropriate time delays between the NE 
and SW lobes with respect to the observer/core 
reference frame can be determined using the 
known lobe velocity j3 and orientation 9 to the 
line of sight. This model predicts the delay of 
the lobe flares tme and t$w, with respect to the 
driving core flare as viewed by the observer, to be: 

TNE - (tl - t jP ( / g._ / g) ' (Ld > 
N „(1 + 0nCOS6) 

tsw = (ti - to)/? ■ (14) 

Here to is the time of ejection of the lobes from the 
core, and t\ is the time of a core flare. With this 
model, the NE lobe should flare at time t\ + tne 
and the SW lobe at time t\ + tsw later than the 
core, as viewed by the observer. The term (ii — 
t ) (3 is the true distance of the lobes from the core 
at time t\. 

Figure 16a shows the observed flux density ver- 
sus time for the three components. Figure 16b 
shows the flux density correlations after correct- 
ing for the model delay, assuming f3j = 1. Figure 
16c shows the same data assuming (3j = 0.90 and 
the correlations between the core with the NE and 
SW lobe flux density flares are significantly worse. 
Additional discussions of the correlations and its 
relationship to various properties of Sco X-l were 
presented in Paper I. 

A lower limit to the flow speed in the beam, 
as indicated by the difference in correlation be- 
tween Figures 16b and 16c, can be determined 
more quantitatively It seems unlikely that the ap- 
parent correlations in Figure 16b are pure chance. 
There are three pairs of flares (C3-L3, C4-L4, C3- 
S3) which are clearly correlated; all are associated 
with a significantly different delay using a model 
with only one free parameter, derived solely from 
the positional information of the components. The 
correlation of the flares for f3j — 1.0 is clearly bet- 
ter than that for fij = 0.9 (and by symmetry for 
(3j = 1.1), but how much better? 

An estimate of the uncertainty can be obtained 



by determining the speed necessary to align each 
of the three flare correlations (C3-N3), (C4-N4) 
and (C3-S3). These results yield the three inde- 
pendent estimates of of 1.10 ±0.10, 1.00 ±0.04, 
and 1.02±0.04 respectively, which gives a weighted 
average /3j = 1.02 ± 0.04. Changing the space 
velocity of the lobes within the estimated uncer- 
tainty affects [3j only at the few percent level. 
The relatively large beam speed of 1.10c associ- 
ated with flare (C3-N3) may be caused by an ad- 
ditional delay in the flare of the core. We have 
suggested that the core may be at the base of the 
jet, at a distance of one or two mas from the bi- 
nary. This would cause a delay in the core flare 
reaching the observer by w 0.1 hours, compared 
the event in the binary system. 

Thus, we derive a lower limit of pj > 0.95 based 
on the above analysis. Although the analysis indi- 
cates that this limit is significant at the two-sigma 
level, the model may be too simplistic. We be- 
lieve that the limit has a confidence level of 70% 
(one-sigma). 

5.9. Flux Density Variations and Doppler- 
Boosting 

Based on the space motion of the lobes, we cal- 
culate that Doppler boosting will produce a flux 
density ratio of the SW lobe to the NE lobe of 
R=0.17 assuming a uniformly filled lobe with an 
isotropic magnetic field. The observed flux den- 
sity ratio is difficult to determine because of the 
variability of the lobes and the time delay in their 
reference frames. However, an estimate can be 
obtained by two methods. First, Figure 13 shows 
the distribution of the flux densities of the NE 
and SW components over all of the observations. 
At 5.0 GHz the highest flux densities observed for 
the NE component are about 18 mJy. For the 
SW component, the largest flux densities are 4.5 
mJy at 1.7 GHz, which corresponds to about 2.2 
mJy at 5.0 GHz assuming a = —0.6. This ra- 
tio is R=0.12. A discussion of the core and its 
relationship to the x-ray and optical emission is 
given elsewhere (Geldzahler et al. 2002). Second, 
another estimate of the Doppler-boosting can be 
made from the flare peak #3 for the NE lobe and 
the SW lobe, shown in Figure 16. The increased 
flux densities at 5 GHz in the flare, compared with 
the baseline derived from a few hours before and 
after the flare, are: Core flare = 3.5 mJy, NE lobe 
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= 17.0 mJy, and SW lobe = 1.0 mJy. The ra- 
tio between the SW and NE lobe flare height is 
R=0.06. This ratio assumes that the energy flow 
in the twin-beams are the same. Finally, hydrody- 
namic modeling of a source with similar properties 
as Sco X-l (Mioduszewski et al. 1997) gives sug- 
gests a ratio of R=0.07 for the peaks of the sim- 
ulated lobes-which is resolution dependent. All 
of these estimates are lower than the R=0.17 ex- 
pected purely from the space motion; however, 
any complexity in the emission properties of the 
lobes (eg. internal structure, non-isotropic mag- 
netic field) tends to decrease this expected ratio 
(Cawthornc 1991) 

5.10. Comparison with other Galactic-Jet 
Sources 

Why does Sco X-l show a lobe phenomenon 
when the observations of other well studied 
Galactic-jet sources generally show multiple com- 
ponents? The general consensus from the other 
well-studied Galactic-jet objects like GRS1915+105 
(Dhawan et al. 2000, Fender et al. 1999), GRSJ1655 
(Tingay et al. 1995; Hjellming & Rupen 1995), 
CygX-3 (Geldzahler et al. 1984; Marti et al. 1992. 
Schalinski et al. 1993; Mioduszewski et al. 2001), 
SS433 (Margon 1984; Spencer 1984), and Cir X-l 
(Fender et al. 1998), is the radiating material is 
probably associated with discrete clouds or shocks 
being energized within a collimated flow of en- 
ergy. Most of these objects are transient X-ray 
sources, preferentially observed just during or af- 
ter a strong flare. These flares may be associ- 
ated with the ejection of a large amount of mate- 
rial that is entrained in a collimated flow of en- 
ergy, producing strong radio emission. However, 
there are several periods when the core of Sco X-l 
showed changing radio emission usually within 5 
mas towards the northeast. Two of these events 
are associated with the emergence of the new NE 
component on MJD 51341.0 (Figure 2) and MJD 
50682.0 (Figure 11), and have already been dis- 
cussed. After the core flare at MJD 51341.95 
(Figure 2), extended emission to the NE persisted 
for about four hours. A compact NE component 
more than 15 mas from the core was also present. 
Figure 14 shows that the core on MJD 50871 was 
also slightly extended by about 3 mas toward the 
north-east direction. The extended emission near 
the core may be associated with material in the 



beam flow. 

The degenerate object associated with Sco X-l 
is probably a neutron star, although the improved 
orbital parameters, based on the orientation of 
the source discussed in §4.2, do not confidently 
constrain the degenerate star mass to less than 
about 2 M Q . For most other Galactic-jet sources, 
the degenerate star is probably too massive to be 
a neutron star. Sco X-l is also a persistent X- 
ray source, unlike many of the other well-studied 
Galactic transient X-ray sources. The relationship 
between the X-ray properties with the nature of 
the degenerate star and with the properties of the 
radio emission is still unknown. 

Finally, the major difference between Sco X-l 
and the other sources may be the viewing angle of 
the source as seen by the observer. The simulated 
images (Mioduszewski et al. 1997), which showed 
a radio structure structure surprisingly similar to 
that seen in Sco X-l, take on the appearance 
of the other Galactic-jet sources when the view- 
ing angle is decreased. At small viewing angles, 
the shocked material moving near c are greatly 
Doppler boosted and dominate over the emission 
from the more slowly-moving lobes. Thus, like 
the extra-galactic sources, the radio appearance of 
Galactic-jet sources may also be a strong function 
of viewing angle. 

6. Comparison with Extragalactic Sources 

One of the major reasons for studying Galactic- 
jet sources is to understand better the jet phe- 
nomenon in extragalactic radio sources. The 
appropriate scaling laws come from considera- 
tions of the accretion disks and infalling mass 
(Begelman et al. 1984). The radio emission is a 
second order product of the accretion process, but 
should also follow the scaling. The luminosity ra- 
tio of Sco X-l to Cygnus A, for example, is 10~ 7 . 
The mass ratio of the neutron star in Sco X-l and 
the massive object in Cygnus A is 10 -9 . We will 
adopt as a typical scaling parameter x, from a 
radio galaxy to Sco X-l x « 10~ 8 . The luminos- 
ity, the time scale and the size will scale as x, the 
density as x^ 1 , the magnetic field as x~ ' 5 , and 
velocity and temperature are invariant to first or- 
der. The following are some examples of simple 
scaling of the properties of Sco X-l's radio lobes 
to those in extragalactic radio sources. 
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Two of the best studied extragalactic ra- 
dio sources with strong lobes are Cygnus A 
(Carilli et al. 1991) and Pictor A (Perley et al. 
1997). Comparisons of these and other sources 
with Sco X-l indicates: 

• Lobe Variability. The lobes in Sco X-l 
vary on an hour time scale. This scales to 
1 x 10 4 yr for radio galaxies. Variability 
in extragalactic lobes have not yet been ob- 
served over about a 30-year span. 

• Lobe Lifetime. The lifetime for a pair of 
lobes in Sco X-l is about one day. This 
scales to 5 x 10 5 yr for a typical extragalactic 
source, and is about the estimated age of the 
Cygnus A hot spots, based on the curvature 
of their radio spectrum (Carilli et al. 1991; 
Cox et al. 1991). 

• Lobe Regeneration. In Sco X-l there 
appears to be replenishment of the ra- 
dio lobes. That is, after a pair of lobes 
fades away, another pair emerges from the 
core, probably in less than one day. The 
one well-documented disappearance of the 
lobes (MJD 51840) occurs after an X-ray 
flare and the emergence of an optically- 
thick core component. The suggestion 
that the the beam flow has been disrupted 
and later reforms is discussed elsewhere 
(Geldzahler et al. 2002). Although a short 
lifetime of 10 5 yr for hot-spots in radio galax- 
ies has been suggested from lobe modeling 
(Cox et al. 1991) and observations of radio 
galaxy lobes (Black et al. 1992), new hot- 
spots are believed to appear in the estab- 
lished radio-lobe region rather than being 
reformed near the core. Hence, unlike Sco 
X-l, the beam flow in extragalactic sources 
may not be completely disrupted. 

• Lobe Size. The linear size of the lobe in 
Sco X-l of 2 x 10 8 km scales to 500 pc for 
radio galaxies, similar to that observed in ra- 
dio galaxies (Black et al. 1992). The west- 
ern lobe in Pictor A (Perley et al. 1997) is 
well-resolved and is somewhat more elon- 
gated perpendicular to the presumed input 
energy flow, as found for Sco X-l. 

• Lobe Distance from Core. The lobe 
separation of Sco X-l is about 40 mas or 



1.6 x 10 km. This scales to a linear size of 
1.6 x 10 18 km or 50 kpc for the typical size of 
a radio galaxy. 

• Spectral Index. The spectral indices asso- 
ciated with the lobes of Sco X-l are similar 
to those seen in extragalactic sources, where 
the most confined lobes have a ss —0.5 and 
the less confined have steeper spectral index. 

• Opening angle. The ratio of the lobe size 
to the distance from the core is probably 
related to the angle opening in the beam. 
For Sco X-l this angle is typically about 1°. 
A similar ratio is observed for extragalac- 
tic double sources that have prominent hot- 
spots in their radio lobes. 

• Magnetic Field Strength. The equipar- 
tition magnetic field in Sco X-l was ~ 1 G. 
If adiabatic expansion of the plasma in the 
lobes is not occurring, the field needed to 
decrease the electron emitting lifetimes from 
synchrotron losses to one hour was ~ 300 G. 
These two magnetic fields strengths scale to 
1 x 10~ 4 G (which is similar to the equipar- 
tition field in Cygnus A) and 3 x 10~ 2 G 
respectively. The configuration of the mag- 
netic field in the lobes of Pictor A appears to 
be circumferential (Perley et al. 1997). The 
degree of polarization in Pictor A is also 
high, indicating an ordered field. No lobe 
polarization information is available for Sco 
X-l. However, if the lobes are dynamically 
controlled by a large magnetic field, a radial 
magnetic field would be more likely with Sco 
X-l. 

A major difference between the radio lobes of 
Sco X-l and those associated with radio galaxies is 
the apparent simplicity of the Sco X-l lobes, usu- 
ally consisting of two hot-spots and no extended 
emission. In contrast, radio galaxies often show 
an extended lobe, with a cocoon of emission, and 
a hot-spot. Emission from the beam between the 
core and the lobe is also observed in some cases. 
Often a secondary hot-spot is observed in a radio 
lobe, not far from the more compact one. These 
secondary hot-spots are thought to be caused by 
either the splatter of material from the primary 
hot-spot, or the remnants of older hot-spots which 
are still radiating (Cox et al. 1991). 
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The difference in appearance between Sco X-l 
and a typical radio galaxy is not the result of the 
high dynamic range images available for many ex- 
tragalactic sources. Even at the modest dynamic 
range of the Sco X-l images (peak flux density 
to rms noise level of about 50) nearly all double 
radio sources have more complicated lobes than 
that observed in Sco X-l. Perhaps the suggested 
strong magnetic field in Sco X-l confines the rel- 
ativistic fluid to a small region. In contrast, on 
MJD 50872 and MJD 51054 the NE lobe of Sco 
X-l resembled the morphology of that in FRII- 
type radio galaxies. The lobe was extended along 
the source axis and the hot-spot was not domi- 
nant. For the MJD 50872 observations, we know 
that the lobe was formed more than 1.5 days ear- 
lier, suggesting that the confinement of relativistic 
particles within a small region is less likely as Sco 
X-l evolves. 

Explanations for the differences seen between 
the radio lobes in Sco X-l and those of extragalac- 
tic sources require a better understanding of the 
formation of the working-surface and the relativis- 
tic particles and fields and their transport into the 
hot-spot. One significant difference between Sco 
X-l and radio galaxies may be the advance speed 
of the hot-spot. For Sco X-l, the advance velocity 
is 0.5c. It is believed that the advance velocity 
in radio galaxies is only mildly relativistic w 0.1c. 
Another difference may be that the lobes in cx- 
tragalactic sources are probably confined by ram 
pressure. In Sco X-l, the lobes may be confined 
by a strong magnetic field, or, not confined at all 
with adiabatic expansion occurring outside of the 
working surface region. 

7. Conclusion 

These extensive VLBI observations of Sco X- 
1 are among the most detailed for a Galactic- 
jet radio source. The high-resolution images over 
the four-year study show relativistic motion and 
variability of all components and a recurrent phe- 
nomenon of lobe production. The resolution and 
sensitivity are sufficient to determine accurate pa- 
rameters of the discrete components as a function 
of time. Since we have observed Sco X-l with- 
out regard to its X-ray and radio flux density, we 
may be sampling the more stationary properties 
of a Galactic-jet source rather than those associ- 



ated with the evolution during the more explosive 
periods. 

The relatively non-speculative conclusions from 
these observations are the following. 

• Sco X-l is composed of three components: a 
radio core near the compact component, a 
NE component which moves away from the 
core, and a weaker SW component which is 
detected about half of the time. 

• Occasionally weak extended emission em- 
anates from the core in the north-eastern di- 
rection, but generally only a NE and SW 
component are detected. 

• All components are variable on time scales 
of about one hour. The rapid increase 
and somewhat less rapid decrease in the 
flux density for all of the components has 
a characteristic time scale of about three 
hours. This time-scale of the radio vari- 
ations in Sco X-l has already been noted 
(Hjellming et al. 1990; Bradshaw et al. 1997) 

• A pair of moving components persists up 
to about two days. Recurrent generation of 
pairs occurs often enough that a naked core 
component is not commonly observed. 

• The comparison of the radio properties be- 
tween the NE and SW components strongly 
suggest that they are intrinsically similar, 
but differ from the effects of relativistic aber- 
ration and Doppler beaming. 

• The average speed of the NE and SW com- 
ponents is 0.45c ± 0.03c, radially away from 
the core at an angle 44° ± 7° to the line of 
sight. 

• The component speeds remain unchanged 
for many hours although the speeds for dif- 
ferent pairs of components at different times 
range between 0.31c and 0.57c. 

• The components of Sco X-l lie on a axis 
which has not varied by more than about 
4° over five years. A variation of axis direc- 
tion of w 3°, with a period of a few days, 
may be present. 
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• The minimum angular size of the NE (and 
presumably the SW) component is 1.5 x 0.7 
mas (6 x 310 1 3 cm) with an orientation per- 
pendicular to the motion of the lobe. This 
corresponds to a volume of 6.1 x 10 40 cm 3 . 
The equipartition magnetic field is between 
0.3 to 1.0 G depending on the ratio of elec- 
trons to protons. 

• The lobes of Sco X-l arc dominated by a 
compact hot-spot. The lobes of most extra- 
galactic sources often contain one or more 
hot-spots and much more extended struc- 
ture. 

More speculative conclusions are: 

• The NE and SW components are lobes gen- 
erated from a working surface where the en- 
ergy flow within a twin-beam from the bi- 
nary impinges on the ISM. This conclusion is 
based on the variability of the components, 
the lack of a systematic decrease in the lobe 
flux density with time, and the simplicity of 
the components. 

• If the above suggestion is correct then the 
temporal correlation of the core flare times 
with those of the NE lobe and SW lobe sug- 
gest a beam velocity > 0.95c. 

• In Sco X-l the lobe emission dominates. 
In most other Galactic-jet sources, emission 
from material or shocks within the jets dom- 
inate. This may be an intrinsic difference 
between the sources, perhaps related to ei- 
ther the nature of the degenerate star or the 
persistence of the X-ray emission. Alterna- 
tively, the difference may be associated with 
the angle between the observer and the di- 
rection of source motion. 

• The variability of the lobes and their linear 
size suggest a radiative lifetime of less than 
one hour. For synchrotron losses in the emit- 
ting plasma, a magnetic field of ~ 300 G is 
needed, which implies a lobe dominated by 
magnetic energy. Alternatively an adiabatic 
expansion of the ultra-relativistic plasma, af- 
ter generation at the working surface, is also 
consistent with the lobe emission character- 
istics. 



• The simple pressure balance in the ultra- 
relativistic plasma on both sides of the shock 
at the working surface predicts an advance 
speed of the lobe in the range of that ob- 
served. The relatively constant speed of the 
lobes over many hours is surprising and no 
explanations are offered. 

• The lobe looks like a wagon wheel. The 
working surface is at the hub and the 
electrons diffuse with 0.57c away rom the 
working-surface shock. The advance of the 
lobe over the lifetime of the radiating elec- 
trons (30 minutes) also contributes to its 
size inthe direction of motion. 

• The slow expansion of the lobes, with a de- 
crease in flux density and steepening spec- 
trum, may be associated with a decreasing 
beam flow and exlarging of the working sur- 
face, as well as from adiabatic expansion of 
material in the lobes. 

• The 'scaling' of the Sco X-l phenomena to 
that of radio galaxies suggests that hot-spot 
variability and life-times may be only 10 5 yr 
in radio galaxy lobes. 

These observations show that VLBI resolutions 
are sufficient to determine the internal properties 
of Galactic-jet radio sources. The time scales of 
the evolution vary from less than one hour to a 
few days or months in some cases. Hence, de- 
tailed, long-term, multi-frequency and reasonably 
continuous observations are required to follow and 
understand the emission properties of the sources. 
Because of the possible dominant effect of the 
magnetic field, linear polarization imaging would 
be very useful. 
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Table 1 

The VLBI Observations of Sco X-l 



Date 


MJD 


UT 


range 


Array 


Freq 


BW 


Flux Dens. 














(GHz) 


(MHz) 


(mjy) 


19-Aug-1995 


49948 


23.7 


to 


26.7 


VLBA 


4.98 


64 


0.4 


16-Mar-1996 


50158 


10.1 


to 


13.0 


VLBA 


4.98 


64 


0.6 


14-Sep-1996 


50340 


22.1 


to 


25.0 


VLBA 


4.98 


64 


4.5 


03-Aug-1997 


50663 


23.9 


to 


29.3 


VLBA 


4.98 


64 


3.5 


21-Aug-1997 


50681 


22.8 


to 


27.8 


VLBA+Y1 


4.98,1-67 


128 


15.0 


27-Feb-1998 


50871 


10.3 


to 


16.0 


VLBA+Y1 


4.98,1.67 


128 


22.0 


28-Feb-1998 


50872 


10.2 


to 


16.0 


VLBA+Y1 


4.98,1.67 


128 


5.0 


29-Aug-1998 


51054 


23.8 


to 


27.0 


VLBA+Y1 


4.98,1.67 


128 


9.0 


ll-Jun-1999 


51340 


02.2 


to 


10.2 


VLBA+Y27 


4.98,1.67 


128 


20.0 


ll-Jun-1999 


51340 


10.2 


to 


18.2 


APT 1 


4.98 


64 


8.0 


ll-Jun-1999 


51340 


18.2 


to 


26.2 


EVN+ 2 


4.98 


128 


18.0 


12-Jun-1999 


51341 


02.2 


to 


10.2 


VLBA+Y27 


4.98,1.67 


128 


14.0 


12-Jun-1999 


51341 


10.2 


to 


18.2 


APT 1 


4.98 


64 


12.0 


12-Jun-1999 


51341 


18.2 


to 


26.2 


EVN+ 2 


4.98 


128 


14.0 


13-Jun-1999 


51342 


02.2 


to 


10.2 


VLBA+Y27 


4.98,1.67 


128 


8.0 



Note. — 

Yl = One VLA telescope; Y27 = The entire VLA 

The APT consists of Australian Telescope Compact Array, Ceduna, Hartebcesthoek Radio 
Astronomical Observatory, Hobart, Kashima, Parkes, Shanghai. 

2 The EVN+ consists of Effclsberg, Green Bank 140', Hartebeesthoek Radio Astronomy 
Observatory, Jodrcll Bank, Mcdicina, Noto, Onsala, Westerbork. 
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Fig. 1. Images at 1.7 and 5.0 GHz with 10 x 5 mas Resolution during 1999 June Observation. 

Five images, rotated 36° counter-clockwise in the sky, at 1.7 GHz (left) and 5.0 GHz (right), are shown for 
each seven-hour period on three consecutive days. The time of each snapshot is shown between the columns. 
An unresolved component will have a position angle of 36°. The minimum contour level (L) and peak flux 
density (P) in mjy are shown next to each plot. The contour levels are — 1, 1, 2, 4, 8 . . . times the minimum 
level. The tick mark are separated by 10 mas on the abscissa and 5 mas on the ordinate. The scale is slightly 
different on the two axes. The dashed lines show the location of the core, NE and SW components. 
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MJD 51342 




-0U6d 
P-J.L 

-0 06 4d 
L-LO 

n njd 
r^L.L 

P*L.] 

r^].L 
F-J.L 
P^J.L 
F-J.J 

L43 

41 Jd 
F-J.] 
L-L ] 




ft 

i. 



J — i 1 - 



ft 



_LLa L| I 



I 



_l I I u_ 



_l I I U I 



_l I I U_ 



I 

-I— 



_1 | | u 



Fig. 2. — Images at 5.0 GHz with 4x2 mas Resolution during 1999 June Observations. Thirteen 
images of Sco X-l, rotated 36° clockwise, are shown for each of the three observing days. The time of each 
snapshot, the peak flux density (P) in mJy and the minimum contour level (L) are given next to each plot. 
The contour levels are -1, 1, 2, 4, 8 ... times the minimum level. The tick marks are separated by 5 mJy, 
except for the abscissa on MJD 51342 where they are separated by 10 mJy. The dashed lines show the 
location of the core, the NE component and the SW component on MJD 51342. The images from the EVN 
and APT arc labeled. Others come from the VLBA. 
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Fig. 3. — Flux Density During 1999 June Observations, (a) Total Flux Density from VLA, (b) Total 
Flux Density from VLBI, (c) Core Flux Density, (d) NE Component Flux Density (e) SW Component Flux 
Density. The flux densities at 1.7 GHz are shown by the (o) points, 5.0 GHz by the (•). The flux density 
scale is logarithmic, hence the flux density difference between frequencies is proportional to the spectral 
index. The one-sigma error bars are indicated. 
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Fig. 4. — The Spectral Index During the 1999 June Observations for the Three Main Compo- 
nents, (a) Core, (b) NE Component and (c) SW component. The spectral indices are calculated only where 
a detection at 1.7 GHz and 5.0 GHz was made. The error bars indicate the one-sigma error. 
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Fig. 5. — Images During 1998 February Observations. The upper part of the figure shows seven 
snapshot images on MJD 50871, rotated by 36° on the sky, at the frequency and resolution indicated. The 
minimum contour levels are 0.5, 0.4, 0.8 mJy for columns 1, 2 and 3 respectively, with contour levels at 
— 1, 1, 2, 4, 8 . . . times the minimum level. The tick marks are separated by 10 mas on the abscissa and 5 mas 
on the ordinate. The lower part of the figure shows an image of Sco X-l on the following day at 1.7 GHz 
and 5.0 GHz, each with a minimum contour level of 0.2 mJy. The image scale is the same as for the previous 
day but the tick marks on the abscissa are 20 mas. 



27 





10 


7 'J a 


i i i i i i i 

■* 






5 


h ® 

% 


(a) Core 


A 

* \ 




1 






f 


>—> 










B 


30 


©a « 
a » s 






>, 


10 




* (b) NE 




Dens 


3 
1 








X 










3 










Oh 


5 




(c) SW 






1 


ftt 
- t| 


i , , , , i 


■>1 ' 



0.5 1.0 1.5 



MJD - 50870 (days) 



Fig. 6. — Flux Density During 1998 February Observations for the Three Components, (a) Core, 
(b) NE Component and (c) SW component. Flux densities at 1.7 GHz are shown by the (o) points, 5.0 GHz 
by the (•) points. The flux density scale is logarithmic and one sigma error bars are indicated. The points 
at each frequency are separated by 50 minutes. 
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Fig. 7.— The Velocity of the NE and SW Components, (top) 1999 June observations and (bottom) 
1998 February observation: The points show the separation of the NE and SW components from the core 
component during the observations. The best linear velocity fit for NE component #1 and component #2 
in 1999 June and for the NE component in 1998 February are indicated by the solid lines. The dashed lines 
through the SW components are exactly 50% of that of the NE component at the same time. 
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Fig. 8. — The Position Angle of the NE and SW Components, (a) The June 1999 observations and 
(b) the Feb 1998 observations. The position angle and error estimate of separation from the core for the NE 
components (•) and SW components (o) are plotted separately. Position angles from both the 1.7 GHz and 
the 5.0 GHz observations are plotted. The dashed-line is at the weighted average position angle 54°. 
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MJD = 49948 MJD=50158 




Fig. 9.— Images at 5.0 GHz for the first four observation, (a) MJD 49948, (b) MJD 50158, (c) 
MJD 50340 and (d) MJD 50663. The peak flux density and lowest contour level are shown in each diagram, 
with contour levels at -1,1,2,3,4,6,8,12 times the lowest level. The cross shows the expected location of the 
core position as determined from the radio parallax and proper motion analysis from all observations. The 
diagonal line is at position angle 54° . Identification of the NE and SW components (sometimes tentative) 
are indicated. The resolution of each image is 6 x 2 mas. 
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Fig. 10. — Images on August 29, 1998. (a) 4.5 x 1.5 mas resolution image at 5.0 GHz, (b) 10 x 5 mas 
resolution image at 5.0 GHz and (c) 10 x 5 mas resolution image at 1.7 GHz. The peak flux density and 
lowest contour level are shown in each diagram, with each diagram indicated the resolution. 
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Fig. 11. — Images on August 22, 1997. (left) 9x3 mas resolution images at 1.7 GHz, (middle) 9x3 
mas resolution images at 5.0 GHz and (right) 3x1 mas resolution images at 5.0 GHz. All images have 
been rotated 36° counter-clockwise. The contour levels are — 1, 1, 2, 4, 8 . . . times the minimum contour level 
of 0.5, 0.4, 0.7 mJy from left to right respectively. The tick mark separation is 5 mas in each coordinate. 
The vertical dashed line on the left of each frame is the core position. The slightly non- vertical dashed line 
indicated the approximate location of the NE component. The dashed line to the left in each frame is the 
location of a possible faint relic NE component. 
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Fig. 12. — The Distance Ratio for the SW and NE Components. The ordinate is the ratio of 
the core-SW component separation divided by the core-NE component separation. The abscissa is the total 
separation between the NE and SW components. The observation date (MJD=T) is indicated for the clumps 
of points, and the error bars indicate the estimated error. The ratio 0.51 ± 0.02 is indicated by the dashcd- 
lines. This value is the mean value and error estimate, heavily weighted by the 1998 February and 1999 June 
Observations. 
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Fig. 13. — The Spectral Index vs Flux Density for the Components, (a) The NE component at 5.0 
GHz, (b) The SW component at 1.7 GHz and (c) The Core at 5.0 GHz. The spectral index a (S oc v a ) is 
plotted versus the flux density (S) of each component. Most of these data are from the 1999 June and 1998 
February observations. The dashed line is at a = —0.6 for (a) and (b), and a = +0.2 for (c). The error bars 
show the estimated spectral index and flux density errors for each point. The data are shown at 1.7 GHz 
for the SW component since it was detected more often at this frequency. 
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Fig. 14. — The Expansion of the NE component on MJD 50871. The position and full- width half- 
intensity angular size of the NE component and core component on MJD 50871 are shown by the ellipses. 
The components have been moved closer together by 22.3 mas in position angle 54° to fit on the plot. The 
numbers at the center of the NE component ellipses indicate the snapshot number (SS) which are separated 
in time by 50-minutes. No snapshot number is given for the core component since it is more stationary. The 
arrow labeled 'c in 50 min' shows the distance that a light signal will travel in the plane of the sky in 50 
minutes. For each snapshot (SS), the total flux density (FLUX) and spectral index (SI) is given for the NE 
component. 
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Fig. 15. — Schematic of Lobe in Sco X-l. A beam of energy from the binary region of Sco X-l interacts 
with the ambient material to form a working surface composed of an ultra-relativistic plasma. Radiating 
electrons diffuse with v = 0.57c along the radial magnetic field lines to form the lobe. The working surface 
is also moving at « 0.45c at 45° to the line of sight. The minimum diameter of the lobe is 1.5 mas (3 x 10 8 
km) and 0.6 mas (1.3 x 10 8 km) perpendicular and parallel, respectively, to the lobe motion. 
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Fig. 16. — The Correlation of Variations in the Core and Lobes in 1999 June, (a) The observed 
flux densities of the core and both lobes, (b) the flux densities with the NE and SW lobes corrected for a 
beam speed of (3j = 1.0 and (c) the flux densities with the NE and SW lobes corrected for a beam speed of 
(3j = 0.9. The flux density scale is logarithmic. The major lobe flares for the core are CI, C2, C3, C4; for 
the NE lobe are Nl, N2, N3, N4; for the SW lobe, the only detected flare, is indicated by S3. The arrows 
on the abscissa of (a) indicated the time of minimum light for the binary system. 
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